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Abstract 

A new strategy for modeling the land surface component of the climate 
system is described. The strategy is motivated by an arguable deficiency in 
most state-of-the-art land surface models (LSMs), namely the disproportion- 
ately higher emphasis given to the formulation of one-dimensional, vertical 
physics relative to the treatment of horizontal heterogeneity in surface prop- 
erties — particularly subgrid soil moisture variability and its effects on runoff 
generation. The new strategy calls for the partitioning of the continental sur- 
face into a mosaic of hydrologic catchments, delineated through analysis of 
high-resolution surface elevation data. The effective “grid” used for the land 
surface is therefore not specified by the overlying atmospheric grid. Within 
each catchment, the variability of soil moisture is related to characteristics 
of the topography and to three bulk soil moisture variables through a well- 
established model of catchment processes. This modeled variability allows 
the partitioning of the catchment into several areas representing distinct hy- 
drological regimes, wherein distinct (regime-specific) evaporation and runoff 
parameterizations are applied. Care is taken to ensure that the deficien- 
cies of the catchment model in regions of little to moderate topography are 
minimized. 
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1 Introduction 

1.1 History of LSM Development 

Atmospheric general circulation models (GCMs) are invaluable tools for ex- 
amining the mechanisms behind climate variability and for characterizing the 
sensitivity of climate to anthropogenic forcing. Of course, to be of value, a 
GCM must produce a realistic climate, and this requires realistic" represen- 
tations of the relevant physical processes. Over continents in particular, a 
proper representation of the land surface energy and water balance is critical. 
The land surface model (LSM) used with the GCM must, for example, prop- 
erly partition the net incoming radiative energy into latent heat flux, sensible 
heat flux, ground heating, and snowmelt energy, and it must properly par- 
tition the precipitation at the surface into evaporation, runoff, and moisture 
storage. Complicating the development of an LSM is the inability, given cur- 
rent computer resources, to resolve explicitly the small-scale structure of the 
physical features that control these partitionings. 

LSMs over the years have varied greatly in structure and complexity. 
In the simplest (and earliest) LSMs, land surface moisture conditions are 
prescribed — dry conditions are artificially maintained in deserts, for exam- 
ple, and wet conditions are maintained in tropical forests. This prescription 
has the advantage of ensuring reasonable evaporation fluxes across the globe, 
but it has the distinct disadvantage of precluding important land-atmosphere 
feedbacks, some of which lie at the heart of continental climate sensitivity. 
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This shortcoming is overcome by “interactive” land surface models. The 
simplest is the so-called “bucket” model of Manabe et al. [1969], which allows 
the water level in a soil moisture reservoir to increase during precipitation 
events and to decrease as the water evaporates. Because the evaporation 
efficiency varies with the water level in the reservoir, rainy periods lead to 
high evaporation rates, and droughts lead to low rates. Runoff occurs only 
when the soil moisture content exceeds field capacity. Manabe’s original 
model has been adapted over the years to include, for example, multiple soil 
layers ( e.g ., Hansen et al. [1983]) and runoff generation from drier soils ( e.g 
Gates and Schlesinger , [1977]). 

The “SVAT” (soil- vegetation-atmosphere transfer) model, a fundamen- 
tally different type of LSM, was introduced in the mid-1980’s by Sellers et 
al. [1986] and Dickinson et al. [1986]. The SVAT approach gives vegeta- 
tion a more direct role in determining the surface energy and water balance, 
particularly by allowing stomatal conductance to decrease in response to in- 
creased environmental stress. The evaporation formulation in most SVAT 
schemes can be described with the aid of a resistance diagram, and in the 
simpler SVAT schemes, this resistance diagram is equivalent to that of the 
Penman-Monteith evaporation formulation [ Monteith , 1965]. The effects of 
vegetation on momentum transfer and on the surface radiation balance are 
also parameterized. 

The roster of models used by participants in a recent major international 
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intercomparison project (the Project for the Intercomparison of Landsurface 
Parameterizations Schemes, or PILPS Henderson Sellers et al. [1993]) indi- 
cates that the basic SVAT approach is the most popular and can, at least 

in that one sense, be considered the “state of the art”. Numerous SVAT 

\ 

variations have been devised; for example, one recent PILPS experiment 
[ Chen et al., 1997] featured about twenty different SVAT LSMs. With some 
modeling groups, the SVAT approach is being enhanced by the inclusion 
of photosynthesis-transpiration physics [Sellers et al, 1996; Dickinson et al, 
1998] and by a treatment of dynamic vegetation [Foley et al, 1996; Dickinson 
et al, 1998]. 

Of particular relevance for this paper are a few recent approaches that 
include an explicit representation of subgrid heterogeneity in surface charac- 
teristics. “Statistical-dynamical” approaches relate the surface energy and 
water fluxes to the spatial heterogeneity of the land surface and the mete- 
orological forcing, characterized statistically [e.g., Entekhabi and Eagleson, 
1989; Famiglietti and Wood, 1990; Avissar, 1992]. The “mosaic” approach 
involves separating the land surface grid element into subgrid tiles based on 
surface characteristics and then treating each tile separately in the simula- 
tion [e.g., Avissar and Pielke [1989], Roster and Suarez [1992], Ducoudre 
et al., [1993], Seth et al. [1994]). A few LSMs have been built around re- 
lationships between topography and the subgrid variability of hydrological 
processes [Famiglietti and Wood, 1991, 1994; Quinn et al, 1995a; Stieglitz et 
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al., 1997; Liang et al., 1994]. As explained in the next section, we consider 
such a focus on subgrid variability to be critical to the realistic operation of 
an LSM. 

1.2 Motivation for a New Strategy 

PILPS [Henderson-Sellers et al., 1993] was designed to document similarities 
and differences in the behaviors of the various state-of-the-art land surface 
models (LSMs) used with numerical atmospheric models. The main finding 
from PILPS was a wide disparity in LSM behavior in a number of environ- 
ments [Pitman et al., 1993; Chen et al, 1997; Wood et al., 1998]. The LSMs 
strongly disagree on how best to simulate land surface energy and water 
fluxes. 

Koster and Milly [1997] addressed these differences with a focused study 
of LSM evaporation and runoff formulations. They avoided comparing spe- 
cific details of the formulations, since an adequately thorough study of this 
type would be unmanageable. Rather, they focused on comparing the un- 
derlying, effective relationships that control the surface fluxes in the various 
PILPS LSMs. In particular, they examined the effective functional relation- 
ships between root zone soil moisture and both evaporation and runoff in 
each model. They found that the annual water balance is largely controlled 
by the shapes of these functional relationships and, in particular, by how 
these relationships are “positioned” relative to each other when plotted in 
the same figure. The upshot of their analysis is that annual evaporation rates 
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are controlled as much by an LSM’s runoff formulation as by its evaporation 
formulation. If the formulation of runoff in an LSM is poor, the LSM will 
produce unrealistic annual evaporation rates regardless of the quality of the 
evaporation formulation. 

* v 

This result may seem obvious, but in some ways it is at odds with the 
apparent development strategies of many modeling groups. The evapora- 
tion and runoff formulations in LSMs are rarely of similar complexity — 
evaporation is usually given much more attention. Standard SVAT models 
account explicitly for canopy interception and for environmental influences 
(e.g., soil moisture availability, air temperature, vapor pressure deficit) on 
canopy conductance, using functions that can be quite complex. In addi- 
tion, as mentioned earlier, a recent trend in model development is focusing 
on explicit treatments of photosynthic controls over transpiration and car- 
bon uptake, and work is progressing in the dynamic prediction of vegetation 
phenology. The one-dimensional nature of the SVAT modeling approach 
encourages these modeling trends by allowing a detailed description of ver- 
tical canopy structure, a framework that is quite adequate for a first-order 
treatment of the essential biological processes involved in transpiration. 

The one-dimensional framework, however, is not amenable to an adequate 
treatment of runoff generation, since runoff in nature is largely controlled by 
spatial heterogeneity in precipitation and surface conditions. Spatial hetero- 
geneity in soil moisture is especially critical. Overland flow by the Dunne 
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mechanism requires rainfall to impinge on a saturated ground surface, typi- 
cally a small fraction of the land area receiving the rain. Overland flow by 
the Horton mechanism is generated when rainfall rate exceeds the infiltra- 
tion capacity of the soil, and this infiltration capacity varies considerably in 
space with both soil texture and water content. Baseflow into rivers and 
streams reflects the three-dimensional structure of the water table. Given 
this complexity, an LSM that relies solely on a one-dimensional soil column 
for its hydrological calculations cannot hope to capture the main physical 
controls on runoff production. Thus, typical LSM runoff parameterizations, 
though given some attention over the years [ e.g Entekhabi and Eagleson, 
1989; Ducharne et al., 1998], are arguably crude in comparison to typical 
evaporation formulations. 

This deficiency suggests a logical development path, namely, an improved 
treatment of the subgrid horizontal structure of land surface hydrological pro- 
cesses. We follow this path here. In our new LSM, subgrid heterogeneity in 
surface moisture state is treated statistically, since computational constraints 
(now and in the foreseeable future) prevent its explicit resolution. Neverthe- 
less, the applied distributions are related sensibly to the topography, which 
exerts a major control over much of the subgrid variability. This represen- 
tation of subgrid soil moisture variability allows, in principle, the explicit 
modeling of the runoff mechanisms outlined above. As an added benefit, it 
improves the representation of evaporation. 
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We describe our new LSM in sections 2 and 3 below; section 2 provides 
a broad overview of the modeling strategy, and section 3 provides specific 
details of the formulation. A discussion of the advantages of this modeling 
strategy follows in section 4. 

The LSM has been fully coded and has been tested against observations. 
Details, along with a discussion of some unique parameter estimation issues, 
are provided in a companion paper [Duchame et al, this issue], hereafter 
referred to as Part 2. 

2 The Catchment Model: Overall Framework 

We now present a nontraditional land surface modeling framework that in- 
cludes an explicit treatment of subgrid soil moisture variability and its effect 
on runoff and evaporation. Many of the ideas presented here and in the next 
section, which describes specific features of the model, are culled from works 
in the literature that have already recognized the importance of these issues. 
These works include those of Famiglietti and Wood [1991, 1994], who pio- 
neered the idea of adding an energy balance to the TOPMODEL formalism 
[Beven and Kirkby , 1979; Beven et al 1994; Sivapalan et al., 1987] and 
applying the resulting LSM to each grid cell of a GCM, Stieglitz et al. [1997], 
who described a computationally efficient method of applying TOPMODEL 
equations to the land surface component of a GCM, and Liang et al. [1994], 
whose VIC model accounts explicitly for subgrid variability in infiltration ca- 
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pacity. The macroscale TOPLATS model of Famiglietti and Wood [1994] is 
indeed similar in concept to the present model in many ways. The two mod- 
els differ significantly in many respects, though, such as in the interpretation 

of soil moisture distribution functions (see section 3.4) and in the variables 

\ 

used to describe soil moisture content above the water table. 

Our description of the model’s overall strategy divides itself naturally 
into two parts: (?) a description of the land surface element used and its 
connection to the overlying GCM atmospheric grid, and (n) a description of 
the treatment of soil moisture variability within the element. 

2.1 Partitioning of Continental Surfaces into Catch- 
ments 

A key innovation in our approach involves the shape of the land surface 
element. We abandon the traditional approach of defining quasi-rectangular 
land surface elements with boundaries defined by the overlying atmospheric 
grid. Instead, we define the fundamental land surface element to be the 
hydrological catchment, with boundaries defined by topography. 

Catchment boundaries have already been established over North Amer- 
ica through the application of a geographical information system to a 30-arc- 
second resolution (approximately 1 kilometer) digital elevation model (DEM) 
provided by the U.S. Geological Survey. The catchment delineation in south- 
western North America, for example, is shown in Figure 1. The delineation 
procedure considered network topology and drainage area [ Verdin and Jen- 
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son, 1996; Verdin and Verdin , 1999] along with the application of ordering 
rules associated with a catchment coding system. 

The average size of the catchments in North America is 3640 km 2 . Thus, 
for a GCM run at a 4°x5° resolution, which in midlatitudes is roughly equiv- 
alent to a grid cell of size 200000 km 2 , the breakdown of the continent into 
catchments represents a significant degree of added resolution at the surface. 
Two catchments lying below the same atmospheric grid cell can be distinct 
in terms of the topographical parameters that describe soil moisture distribu- 
tions (see below) and the vegetation type assumed for the catchment, which 
can be extracted from high resolution vegetation data sets [DeFries, 1994]. 

The mismatch between the regular atmospheric grid and the irregular 
catchment grid makes necessary an algorithm that disaggregates the atmo- 
spheric forcing to the catchment scale and aggregates the catchment “prod- 
ucts” (i.e., the surface turbulent, radiative, and momentum fluxes) to the 
GCM grid scale. For the aggregation of surface or near-surface fluxes onto 
the atmospheric grid, we use a conservative areal weighting scheme. The net 
flux, Fi, into atmospheric grid cell i is computed with 


Fi = 


nF ' n-<4 n j 


( 1 ) 


where A n j is the area of catchment n that lies within the area defined by grid 
cell i. For the transformation of radiative fluxes and precipitation defined 
on the atmospheric grid onto the catchment grid, we use (in the present 
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incarnation of the model) the analogous conservative scheme: 


F n 


FiiFjAni 

^i-A n i 


(2) 


We currently apply (2) to the atmospheric state variables as well. 

Of course, the coupling between the land and atmospheric models can 
(and should) be made more complex than this. Sub-gridscale variations in 
surface temperature and turbulent fluxes can, in principle, have important 
effects on the development of the boundary layer and on the generation of 
mesoscale circulations. Such inferred patterns, and variations in the topog- 
raphy itself, can in turn be incorporated into the disaggregation of the GCM 
forcing. Various disaggregation methods [ e.g Gao and Sorooshian, 1994; 
Leung and Ghan, 1995] are available in the literature and will eventually be 
incorporated into our overall approach. A full discussion of coupling issues 
is reserved for a future paper. 


2.2 Subgrid Variability Within Each Catchment 

As mentioned above, the separation of the continental surface into catchment 
elements allows for an explicit treatment of subgrid-scale heterogeneity at the 
land surface. Our strategy goes much further in this regard, though, because 
soil moisture is assumed to vary significantly within each catchment element. 
This is made straightforward by the use of a pre-existing, well-tested model 
of catchment processes, namely, TOPMODEL of Beven and Kirkby [1979]. 
We use this model to diagnose root zone soil moisture distributions from 
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the morphology of the catchment and from our bulk soil moisture prognos- 
tic variables. (Note that alternative catchment hydrology models are also 
available, and some may be superior to TOPMODEL for describing many 
situations. As our catchment model evolves, we may choose to replace the 
TOPMODEL component with one of these alternatives. In other words, 
our overall framework is flexible; the use of TOPMODEL is simply consid- 
ered a suitable starting point toward a proper representation of subgrid soil 
moisture variability. This is discussed further in section 4.1.) 

The derived distribution of soil moisture in the root zone allows the sep- 
aration of the catchment into specific hydrological regimes. Specifically, we 
derive: (z) the fraction of the catchment over which the ground surface is com- 
pletely saturated (e.g., along riverbeds), (it) the fraction of the catchment 
over which the ground surface is not saturated but transpiration nevertheless 
proceeds at the unstressed rate, and (in) the fraction of the catchment over 
which the soil is too dry to allow any transpiration (i.e., the soil moisture is 
at or below the wilting point). This separation is the catchment strategy’s 
most important advantage over the traditional one-dimensional SVAT model. 
The physical processes controlling runoff and evaporation are very different 
in these different regimes. Thus, the explicit partitioning of the catchment 
into the different areas and the application of different parameterizations in 
each area should lead to more credible estimates of evaporation and runoff 
across the catchment. 
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For example, the catchment LSM shuts off transpiration in the wilting 
fraction and computes unstressed transpiration in the other two fractions. 
Bare soil evaporation naturally proceeds at a very high rate in the saturated 

fraction, at a moderate rate in the unstressed fraction, and at a low rate in the 

\ 

wilting fraction. Infiltration and soil moisture diffusion effectively proceed 
at different rates in the different areas. Also, as in Stieglitz et al. [1997], 
we convert all precipitation into overland runoff in the saturated fraction 
of the catchment and allow part of the precipitation to infiltrate the soil in 
the other fractions. The water table distribution associated with the soil 
moisture distribution allows a physically-based calculation of baseflow. 

At each time step, the surface water and energy fluxes computed for each 
areal fraction are combined into a single catchment flux, and the fluxes from 
adjoining catchments are then combined via (1) for input into the GCM. The 
areal fractions within a catchment, of course, respond dynamically to changes 
in the bulk water variables. Thus, the saturated fraction can grow and the 
wilting fraction might disappear during pluvial periods, and the reverse may 
occur during extended dry periods. The strength of these variations depends 
in large part on the topography. Given that we allow such dynamic changes 
in areal fraction, care is taken to ensure energy and water conservation across 
the catchment at all times. 
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3 The Catchment Model: Specific Details 
3.1 The TOPMODEL Framework 

Again, in its first incarnation, our LSM uses TOPMODEL equations to relate 
the water table distribution to the topography. (TOPMODEL in this paper 
refers to the original formulation of Beven and Kirkby [1979] and not to the 
fully developed LSM of Beven et al. [1994].) The TOPMODEL formulation 
[Beven and Kirkby , 1979; Beven, 1986a, b; Sivapalan et al., 1987] allows a 
dynamically consistent calculation of the spatial distribution of water table 
depth in a catchment from knowledge of topography statistics. At the heart 
of TOPMODEL are three basic assumptions: (i) the water table is nearly 
parallel to the soil surface so that the local hydraulic gradient is close to tan 
(3, where p is the local slope angle; (ii) the saturated hydraulic conductivity 
falls off exponentially with depth; and (Hi) the water table is recharged at a 
spatially uniform, steady rate that is slow enough (relative to the response 
time of the catchment) to allow the assumption of a water table distribution 
that is always at equilibrium. Given these assumptions, an analytic relation 
can be derived between the catchment mean water table depth, d, and the 
local water table depth, d, at an interior point: 

rf=3- H h w' :f )’ (3) 

where In (a/tan/?) is the “topographic index” at the point in question and x 
is the mean catchment value of the topographic index. The term a is the 
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upstream area that contributes flow through a unit contour positioned at 
the point, and u is the parameter describing the vertical profile of saturated 
hydraulic conductivity. 

This framework by itself provides a method for calculating, from a single 

\ 

bulk moisture variable, the saturated areal fraction of the catchment (section 
3.2) and the baseflow from the catchment. It cannot provide everything, 
however. In particular, because it focuses on equilibrium conditions for the 
water table, it does not address non-equilibrium moisture conditions in the 
vadose zone. Such non-equilibrium conditions are common and are critically 
important in the root zone, for example. Thus, the basic TOPMODEL 
framework is amended, as described below. 

The distribution of topographic indices within a catchment can be ex- 
tracted directly from digital elevation models (DEMs). Proper consideration 
must be given to the effects of DEM resolution on the character of the re- 
sulting distribution [Wallace and McCabe 1995; Quinn et al., 1995b; Wolock 
and Price, 1994; Zhang and Montgomery, 1994], as discussed further in Part 
2 [Duchame et al, this issue]. 

3.2 Soil Moisture Prognostic Variables 

3.2.1 Catchment Deficit 

The first bulk variable in our TOPMODEL-based catchment model is the 
“catchment deficit”, Md ■ The catchment deficit is defined as the average 
amount of water, per unit area, that would have to be added to bring all of the 
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soil throughout the catchment to saturation, assuming that the unsaturated 
zone is initially in an equilibrium state. 

This is a rather unconventional prognostic variable, and it thus deserves 
additional explanation. Figure 2 outlines the calculation of Mp. The top 
part of the figure shows the equilibrium profile of soil moisture with depth at 
an arbitrary point in the catchment. The part of the profile above the water 
table is determined from the balance between the pressure head gradient, 
which tends to draw moisture up, and gravity. The precise form of the 
equilibrium profile, as derived from the relations of Clapp and Homberger 
[1978], is 



where w is the degree of saturation (or "wetness”) at a height z above the 
water table, ip s is the matric potential in the soil at saturation, and b is a 
soil parameter. By vertically integrating 1 -w(z) over the unsaturated zone, 
we can determine the local moisture deficit, i.e., the amount of water that 
would need to be added at that point to bring the soil column to complete 
saturation. 

The bottom part of Figure 2 illustrates the intuitive notion that the lo- 
cal moisture deficit varies significantly across the catchment, with low values 
where the water table is near the surface (at the bottom of hills) and high val- 
ues where the water table is deeper (at the top of hills). Given an analytical 
distribution of the topographic index, the TOPMODEL equations provide an 


17 


analytical distribution of water table depth in the catchment, from which an 
analytical distribution of local deficit can be readily derived. The catchment 
deficit can then be calculated as the integral of this local deficit across the 
catchment’s area. 

\ 

Note that this is a somewhat unusual (though not novel) application of the 
TOPMODEL framework. The TOPMODEL equations are not designed to 
generate soil moisture contents in the unsaturated zone and are not derived, 
in any case, from an assumption of equilibrium conditions above the water 
table. We simply use the TOPMODEL equations to describe the shape of the 
water table at any given time and then apply the equilibrium conditions as 
a straightforward way of characterizing, to first order, the moisture content 
above it. The combination of the TOPMODEL water table assumptions 
with the assumption of equilibrium moisture profiles allows a given value of 
catchment deficit to be associated with a unique description of equilibrium 
horizontal soil moisture variability. 

To eliminate an unrealistic drift in the moisture state variables during 
extended dry periods, the model imposes a maximum catchment deficit based 
on estimated soil profile depths [ e.g Webb et al., 1993]. The implications 
of this maximum value on the calculation of subgrid moisture regimes is 
discussed in section 3.4. 
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3.2.2 Root Zone Excess 


If equilibrium conditions could indeed be assumed in the unsaturated zone, 
the catchment deficit by itself would be sufficient to characterize the catch- 
ment’s complete moisture state — no other moisture variable would be 

\ 

needed. Equilibrium conditions, however, are generally not the rule. Soil 
moisture near the surface ( e.g ., in the root zone) responds relatively quickly 
both to the infiltration of precipitation water and to the extraction of water 
via evapotranspiration. Given the critical role played by root zone moisture 
in the surface energy and water balances, we require a moisture variable that 
accounts for such non-equilibrium behavior. For this reason, we introduce 
the “root zone excess”. 

The root zone excess at a given point is defined as the amount by which 
the moisture in the root zone there exceeds (or is less than) the moisture 
content implied by the local equilibrium moisture profile. A positive point 
root zone excess is indicated by the lightly shaded area in Figure 3. The 
value becomes negative, of course, if evaporation exceeds precipitation over 
an extended period of time. 

The point root zone excesses need not be the same throughout the catch- 
ment. The positive excesses that tend to follow a large storm, for example, 
cannot be realized in the areal fraction of the catchment that, at equilibrium, 
is already saturated. Nevertheless, our formulation uses a single, “bulk” root 
zone excess variable that represents the areal integral of the point excesses 
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across the catchment. In other words, the root zone excess for the catchment, 
Mrzi is the average amount of water, per unit area, by which conditions in 
the root zone across the catchment are out of equilibrium. The subcatch- 
ment distribution of point root zone excess is diagnosed from the bulk value 
through considerations of topography, vegetation, and other catchment fea- 
tures. 

Combining the diagnosed distribution of point root zone excess with the 
distribution of equilibrium root zone moisture (as diagnosed from the catch- 
ment deficit) produces a distribution of total root zone moisture. This dis- 
tribution is used to partition the catchment into subregions, as described in 
section 3.4. 

3.2.3 Surface Excess 

The third soil moisture prognostic variable is akin to the “surface soil layer” 
employed by many one-dimensional SVAT models. The “surface excess”, 
•^se> essentially describes the degree to which the water in the top several 
centimeters of soil is out of equilbrium with the water in the root zone below 
it. A given catchment deficit and root zone excess uniquely define a distribu- 
tion of root zone soil moisture, which in turn, through application of (4), can 
be used to diagnose a corresponding equilibrium distribution of surface mois- 
ture. The surface excess is defined as the average amount of water (per unit 
area) by which the actual surface moisture exceeds this equilibrium amount. 

A positive surface excess is induced by the infiltration of rainwater or 
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snowmelt, and a negative surface excess is induced by evaporation from the 
soil surface. The presence of this surface moisture variable was found to be 
crucial to an adequate treatment of bare soil evaporation and surface runoff 

generation. 

\ 

Because the surface excess is assumed not to vary in space, the value for 
the catchment as a whole (i.e., the actual prognostic variable) is the same 
as that assumed for any local point. The total surface moisture nevertheless 
varies in space due to the topographically-induced variations in equilibrium 
soil moisture — the surface soil at the top of the hillslope will be drier than 
that on the valley floor. If the assumption of uniformity leads to any local 
impossibilities in moisture content, such as moistures in excess of saturation, 
the catchment value of surface excess is adjusted accordingly. 

3.3 Transfers Between Moisture Reservoirs 

3.3.1 Transfer Between the Root Zone and the Water Table 

The root zone excess continually interacts with the equilibrium water state. 
If the root zone excess is positive, some water is transferred from the root 
zone to the equilibrium state, and the water table rises. Alternatively, if the 
excess is negative, the root zone gains some water and the water table falls. 
Thus, this interaction, applied at each simulation timestep, always acts to 
reduce the absolute magnitude of the excess — it always acts to bring the 
total system closer to equilibrium conditions. 

The calculation of this moisture transfer, AM, takes advantage of the 
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subcatchment soil moisture distributions diagnosed from the two bulk mois- 
ture variables and the topographic characteristics. To see how this is done, 
consider first a somewhat impractical approach. The catchment is separated 
into a multitude of small pixels. At each pixel, the local root zone excess and 
depth to the water table is established from the subcatchment distributions, 
and the Richards equation is solved (over a highly resolved discretization of 
the vertical moisture profile, using hydraulic properties consistent with the 
basic TOPMODEL formulation) to approximate the local transfer between 
the excess and the underlying equilibrium profile. The pixel-based transfers 
are then summed across the catchment to produce AM, which is then used 
to update the two bulk variables. 

The impracticality of this distributed approach stems from its high com- 
putational expense. The conceptual advantages of a distributed calculation, 
however, can nevertheless be captured by a much simpler formulation. We 
perform, prior to any climate simulation, detailed distributed calculations in 
a given catchment over a wide range of M D and M RZ values. The results are 
processed into a simple, approximate catchment-specific relationship between 
AM, Mp, and M RZ , having the following form: 

A M„ = —M t2 — (5) 

Ti 

and 

A M d = -M„— , (6) 

T 1 
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where /\t is the the timestep length and T\ is a timescale of moisture transfer, 
an empirical function of M& and A/rz- The precise form of this empirical 
function, as derived from the distributed calculations, is described in Part 

2 [ Duchame et al, this issue]. In essence, Ti decreases with decreasing 

■ \ 

and with increasing Mr.z- Its sensitivity to these bulk variables varies with 
topography and thus varies from catchment to catchment. 

3.3.2 Transfer Between Surface Excess and Root Zone Excess 

The magnitude of the surface excess, M se , is reduced at every timestep 
through the transfer of moisture between the surface excess and root zone 
excess variables: 

AM se = — M se — , (7) 

r 2 

and 

At 

A M„ = M se —, (8) 

where t 2 is a timescale analogous to r x in (5) and (6). The value of r 2 
decreases with an increase in either the total root zone moisture or the surface 
excess. The form of the empirical function used for r 2 is presented in Part 2 
[Duchame et al., this issue]. 

3.4 Spatial Partitioning Within the Catchment 

The bulk moisture variables allow the catchment to be separated into three 
distinct regions. The “saturated region” consists of all points for which the 
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root zone is fully saturated, and it has a fractional area of A, at . The “transpi- 
ration” region, consisting of all points having subsaturated root zone mois- 
tures that lie above the vegetation-specific wilting point (0 w iit)i has a frac- 
tional area of /hi- Finally, the “wilting” region, consisting of all points for 
which transpiration is shut off completely, has a fractional area of A w n t . 

The areal fractions are determined through manipulation of a derived 
probability density function (pdf) of root zone soil moisture. If the water 
table lies above the assumed bedrock depth, then associated with each point 
in the catchment is a specific equilibrium moisture profile, represented by 
(4). Integrating this vertical profile across the root zone alone gives a local 
equilibrium root zone moisture, m IZ : 

1 r 5 

m TZ — — _ w(z)dz, (9) 

d TZ Jd-drz 

where d rz is the depth of the root zone. Equations (4) and (9) allow us 
to transform the pdf of water table depth into a pdf of equilibrium root 
zone moisture across the catchment. This pdf has a complex form that is 
nevertheless amenable to a simple approximate representation, as described 
in Part 2 [Ducharne et al . , this issue]. A typical shape of this pdf is shown 
in Figure 4 by curve (A); the pdf is characterized in part by a minimum soil 
moisture, 0 O - 

Again, the pdf described above is for the equilibrium soil wetness. If the 
water table depth corresponding to the current catchment deficit, however, 
lies below the assumed bedrock depth, then no free-standing water table is 
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assumed to exist in the catchment, and the determination of the "equilib- 
rium” pdf (still defined as the pdf associated with the current value of M D ) 
must be revised. This is done by assuming specific shapes for the pdfs at 
three critical values of Mp. These shapes are indicated by the curves (B), 
(C), and (D) in Figure 4. The first pdf, (B), corresponds to the point at which 
the mean water table depth is exactly the assumed depth to bedrock; this is 
taken to be the last point at which an assumed water table can influence the 
shape of the root zone wetness pdf. The second pdf, (C), corresponds to the 
onset of wilting (0„ = 0 w iit)- The shapes of (B) and (C) are the same because 
prior to wilting, the main sink of moisture, evapotranspiration, is assumed 
for convenience to act uniformly on the root zone moisture — unstressed 
transpiration does take place, after all, in both the saturated and transpi- 
ration fractions. Thus, the soil in the non-wilting catchment is assumed to 
dry uniformly, leading to a simple translation of the pdf. (Note that with 
the bedrock depths we assume, which are consistent with the bedrock depths 
typically assumed implicitly by TOPMODEL, the relative positions of (B) 
and (C) are always maintained, i.e., at state (B), the soil moisture everywhere 
lies above the wilting point.) The final pdf we consider, (D), corresponds to 
complete wilting and is represented by a simple delta function centered at 

^wilt • 

The values for ^4 sat , A lT , and j4 w ii t vary according to where Mq and the 
current value of the root zone excess, iV/ rz , position the "non-equilibrium” 
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pdf of soil moisture (i.e., the pdf corrected for root zone excess) in Figure 4. 
M rz basically acts to translate along the 9 axis the pdf established by \I D . 
Details are presented in Appendix A. 

Again, the separation of the catchment into the three subregions consti- 
tutes a critical advantage of the catchment model. The physical mechanisms 
that control evaporation and runoff generation in these subregions are fun- 
damentally different. The explicit separation of the subregions allows us to 
employ the appropriate treatments in each. 

Such manipulation of pdfs, by the way, distinguish this LSM from some 
other topography-based LSMs, such as TOPLATS [Famiglietti and Wood, 
1994] and the model of Quinn et al. [1995b]. In these other models, the 
pdf of topographic index itself is discretized, and calculations are performed 
over each pdf interval. In the catchment LSM described here, the analytical 
topographic index distribution is effectively transformed into an analytical 
distribution of root zone moisture, which in turn is discretized into at most 
three hydrological regimes of relevance. 

3.5 Surface Energy Balance Calculations 

Full details of the energy balance calculations are not provided here, since 
they parallel exactly those in the Mosaic LSM of Roster and Suarez [1992, 
1996]. We focus instead on those aspects of the calculations that are unique 
to the catchment model. In effect, a full one-dimensional SVAT model com- 
putes the energy balance separately in each of the three subregions. In accor- 
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dance with the modeled distribution of soil moisture, the resistances applied 
to the evapotranspiration calculation (a critical component of the energy 
balance calculation) vary significantly between the subregions. For exam- 
ple, resistance to bare soil evaporation, a function of surface soil moisture 
[Koster and Suarez , 1996], is very small in the saturated region, is moder- 
ate in the transpiration region, and is high in the wilting region. The same 
small but nonzero resistance (a function of vegetation type, solar radiation, 
and ambient air temperature) limits transpiration in both the saturated and 
transpiration subregions, whereas in the wilting region, this resistance is set 
high enough to shut down transpiration completely. 

Each subregion maintains its own prognostic surface temperature; no 
“smoothing out” of this temperature is performed at the end of a time step. 
This allows the valley bottoms, where more evapotranspiration occurs, to 
remain consistently cooler than the drier hilltops. 

Note that an intermediate “moisture-stressed” evaporation regime is not 
included in this formulation — transpiration proceeds at an unstressed rate 
in the two wetter subregions and does not occur at all in the driest subregion. 
We neglect moisture-stressed nonzero transpiration for two reasons. First, 
the moisture regime over which it occurs is relatively small; once plants begin 
to be moisture-stressed, only a small additional amount of drying leads to 
total wilting. By ignoring this small regime, we can avoid the application 
of a complex moisture-stressed transpiration model, which requires param- 
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eters that are very difficult to estimate. Second, a key reason for including 
moisture-stressed evaporation in many LSMs is to allow a smooth transition 
between the full transpiration and full wilting regimes. The catchment model 
achieves this smooth transition through the dynamically varying areas of the 
subregions. 

3.6 Ground Thermodynamics 

Although each of the three subregions maintains its own surface/canopy tem- 
perature, temperatures at deeper levels are assumed to be spatially homoge- 
neous. Figure 5 illustrates the framework applied to the calculation of the 
ground temperatures. The net heat flux from layer 1 to layer 2 in the figure is 
computed by weighting the individual heat flux from each surface subregion 
by its respective area and then summing together the three weighted fluxes. 
Heat transport within the soil column is governed by linear diffusion along 
the thermal gradient. As indicated in the figure, layer thicknesses follow a 
geometric series. The top layer’s thickness is taken to be 5 cm to allow us to 
capture the diurnal range in the surface radiating temperature. The ground 
profile extends to a depth of 10m to be compatible with the assumption 
that the bottom boundary of the deepest model layer is completely insulat- 
ing; this depth is approximately three times the seasonal damping depth of 
typical soils. 

Note that the thicknesses associated with the layers in Figure 5 are un- ' 
related to the layer thicknesses associated with the moisture prognostic vari- 
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ables, being chosen instead for consistency with the heat diffusion calcula- 
tions. Even so, heat conductivity does increase in a simple way with increas- 
ing bulk soil moisture. 

3.7 Canopy Interception, Throughfall and Surface Runoff 

\ 

A prognostic moisture variable not mentioned above is canopy interception, 
i.e., the easily evaporated water that sits on leaves and ground litter follow- 
ing rain events. A single interception reservoir overlies all three catchment 
subregions. Interception storage and evaporation is treated exactly as it is in 
the Mosaic LSM [ Koster and Suarez , 1996]. The storage capacity is a simple 
function of leaf area index. A simple representation of subgrid storm po- 
sition “memory” ensures that some of the precipitation (the “throughfall”) 
will almost always fall through the canopy to the ground surface. 

The throughfall, Pj - , falls uniformly on all three catchment subregions. 
Throughfall impinging on the saturated region is immediately converted into 
surface runoff. Throughfall impinging on the other areas infiltrates the soil, 
unless the soil is frozen or the surface excess, M se , is positive, indicating 
a rain-induced build-up of moisture near the surface. In the latter case, 
some overland flow, Q s , is allowed from the transpiration area only. The full 
equation for overland flow is: 

Qs = PrA sat (M^ < 0) 

Qs = Pt (a sat + Ar Se — ) (Mse > 0) (10) 

\ ^se-max / 
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In (10), Mse-max is the maximum possible value of the surface excess given 
the current values of catchment deficit and root zone excess. If infiltration 
saturates the surface layer, the excess over saturation is converted into surface 
runoff. 

Note that (10) is not based on any detailed offline calculation. It is simply 
designed to capture, in a crude way, a smooth transition between zero runoff 
production when M se is zero to a complete conversion of rainfall into runoff 
when the surface is completely saturated. 

3.8 Baseflow 

The calculation of baseflow, or the loss of groundwater to streams, is straight- 
forward given the use of TOPMODEL equations to describe the distribution 
of water table depth. The TOPMODEL framework [Sivapalan et al, 1987] 
relates baseflow, G, directly to the mean water table depth, d: 

G = -^ (sttrface >e-*-3 (11) 

Because the catchment deficit in a given catchment is a unique function of 
the mean water table depth, a simple relationship between baseflow and Mq 
can be easily derived. The particular relationship used is described in Part 
2 [ Duchame et al ., this issue]. 

As discussed in section 3.4, a bedrock depth is assumed in the calculation 
of A sat , A tr , and A wi i t . Baseflow is shut off when d falls below this bedrock 
depth. 
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3.9 Snow Model 


The snow model of Lynch- Steiglitz [1994] has been coupled to the catch- 
ment model. This three-layer snow model accounts for snow melting and 
refreezing, dynamic changes in snow density, snow insulating properties, and 
other physics relevant to the growth and ablation of the snowpack. It has 
been tested with considerable success against detailed measurements in the 
Sleepers River catchment of Vermont. 

The catchment framework, however, necessitated some modifications to 
the original scheme. In particular, we now ensure a smooth transition be- 
tween snow-free and snow-covered conditions in order to capture the gradual 
growth of a snowpack’s spatial extent and to avoid abrupt changes in the 
surface energy balance calculations. The approach used is straightforward. 
We assume a minimum local snow water equivalent of 13mm, a value that 
allows the resolution of the diurnal surface temperature signal yet still pro- 
duces a stable solution with a 20-minute timestep. If a given volume of snow 
falls on a snow-free catchment, that volume is spread uniformly over a frac- 
tion of the catchment so that the local water equivalent at any snow-covered 
point is 13mm. Thus, if the snow falling in a timestep has a total (water- 
equivalent) volume V in m 3 , and if the area of the catchment is A in m 2 , 
then the snow-covered areal fraction, A s /A, is taken to be 


As 

A 


V 

MSA 


( 12 ) 


The snow-covered areal fraction increases similarly as more snow falls until 
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it reaches 1, at which time the local water equivalents across the catchment 
start increasing uniformly. 

Surface energy calculations are performed separately over the snow free 
and snow covered areas. When the fractional coverage is less than one, the 
snow model is represented with a single snow layer, whereas three model lay- 
ers are used when the snow coverage is complete [ Lynch- Stieglitz , 1994], The 
transition between the single layer and three layer representations involves a 
simple conservative redistribution of layer heat and water contents. 

A paper describing this snow model, its interaction with the simulated 
water balance of the catchment, and its validation against available observa- 
tions is currently in preparation (Stieglitz, pers. comm.). 

4 Discussion 

4.1 TOPMODEL Limitations 

Because the catchment model employs TOPMODEL equations, we must 
pay special attention to the behavior of the LSM in regions of low topogra- 
phy, where the advantages of TOPMODEL are lost. In the limiting case of 
perfectly flat terrain, topography-based variations in soil moisture are elim- 
inated. The land surface element is therefore not areally partitioned into 
distinct moisture regimes. Furthermore, each point in the element has the 
same assumed moisture profile in the unsaturated zone and the same as- 
sumed depth to the water table. As a result, the formulation of subsurface 
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moisture transport, which was described in section 3.3 as equivalent to an 
areal integration of distributed Richards equation solutions, reduces in the 
flat case to a single, one-dimensional application of the Richards equation. 
In the flat case, the LSM can therefore be considered equivalent, to a tradi- 
tional three-layer LSM — it uses the one-dimensional Richards equation to 
transport moisture between a surface moisture reservoir, a root zone reser- 
voir, and an “equilibrium” reservoir that mainly reflects moisture conditions 
deeper in the soil. 

A region’s terrain, of course, is never perfectly level. As topographical 
variations decrease, however, the catchment LSM will capture more and more 
of the flavor of a traditional three-layer model. The use of the root zone and 
surface excess variables allows the catchment LSM to function at least as 
well as traditional LSMs in flatter areas, despite TOPMODEL’s limitations 
there. : ; : 

Other limitations of TOPMODEL must also be considered. Key among 
them is the assumption that the water table within a catchment is everywhere 
parallel to the ground surface, which is an oversimplification, especially at 
larger spatial scales. In fact, at such scales (several kilometers or more), sub- 
surface flow may be controlled more by local geology and geomorphology. We 
nevertheless consider the use of TOPMODEL equations to be the appropriate 
first step in the direction away from the traditional one-dimensional “layer” 
approach to modeling subsurface soil moisture transport. It appears, in fact, 
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to be the most sensible step to take at this time given the availability of global 
topographic data, the lack of other relevant data (e.g., bedrock depths) on 
the global scale, and the inability, in any case, to avoid a statistical char- 
acterization of subgrid moisture distributions under current computational 
constraints. TOPMODEL, though perhaps not optimal for the spatial scales 
we are considering, nevertheless captures the critical differences between ups- 
lope and downslope hydrological behavior and should give a useful first order 
description of subgrid soil moisture variability. This is, in fact, supported by 
the applications of the catchment LSM described in Part 2 [ Duc.harne et al, 
this issue]. 

In any case, as mentioned before, the modeling framework is flexible, and 
like all LSMs, the catchment model will continue to evolve. We can certainly 
incorporate more realistic soil moisture characterizations into the model as 
they (and the supporting data) become available. 

4.2 Additional Conceptual Advantages 

A few additional advantages of the catchment approach should be mentioned. 
First, the partitioning of the continental surface into a mosaic of catchment 
“tiles” allows a straightforward calculation of annual streamflow from large 
river basins. The algorithms used to delineate catchment boundaries from 
digital elevation map data [ Verdin and Jenson , 1996; Verdin and Verdin, 
1999] establish flow directions and identify the set of catchments that lie 
within a given larger basin; the resulting annual streamflow for the larger 
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basin can thus be calculated easily from the catchment values [ Ducharne et 

al., 1998]. This is important, of course, for model validation at the annual 

timescale and for the coupling of LSM runoff products to an ocean model. 

The catchment framework also supports the development and use of river 

\ 

routing schemes that emphasize the timing of river water transport, which is 
crucial for model validation at shorter timescales. 

Second, the explicit treatment of subgrid soil moisture distributions is 
amenable to a more tenable validation of modeled soil moisture against ob- 
servations. Such validation with traditional one-dimensional LSMs, though 
given some attention in recent years ( e.g Robock et al. [1995]), is conceptu- 
ally limited by the nonlinear relationship between soil moisture and evapora- 
tion, given a realistic point evaporation function, a single (in the horizontal) 
moisture variable in a traditional layer-type LSM cannot produce a realistic 
average evaporation rate if it must also represent an areally-averaged soil 
moisture content. With the catchment LSM, realistic evaporation rates and 
soil moisture contents can theoretically coexist. In addition to model valida- 
tion, of course, this has important implications for the potential assimilation 
of observed soil moisture data into the model. 

Third, the framework is amenable to the elevation-based disaggregation 
of the grid-scale atmospheric forcing provided by the GCM. We can make 
use of the detailed topographical characterization of each catchment and its 
neighbors to assign, for example, lower near-surface atmospheric tempera- 
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tures to higher elevation catchments or to the higher areas of an individual 
catchment. These elevation-based temperatures could lead to elevation-based 
variations in surface heat fluxes and snowmelt. Variations in shortwave radi- 
ation with topographical aspect could perhaps be included. Procedures for 
disaggregating precipitation by elevation already exist in the literature [ e.g 
Leung and Ghan , 1995]. 

Finally, and perhaps most importantly, the specification of the catchment 
as the fundamental land surface element provides a more direct link to basic 
hydrological science, for which the catchment is a natural unit of study. Our 
model development can take advantage of watershed studies that span several 
decades. The use of the TOPMODEL framework is considered a useful start- 
ing point in this development. Useful features of other catchment hydrology 
models can be incorporated into our new LSM as they are identified. 

5 Summary 

We have described a strategy for improving what appears to be a “weak 
link” in most land surface models (LSMs), namely the characterization of 
subgrid soil moisture variability and its impact on runoff generation. The 
catchment strategy differs from the traditional SVAT strategy in two im- 
portant ways. First, the boundaries of the basic land surface element, the 
hydrological catchment, are derived from topographic data. They are thus 
irregular and independent of the overlying atmospheric grid. Second, sub- 
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catchment distributions of root zone soil moisture are diagnosed from bulk 
moisture variables and topographic characteristics. The distributions are 
used to divide the catchment’s area into three sections, each representing 

a distinct moisture regime. Evaporation and runoff processes are modeled 

* \ 

differently in each section, thereby producing a more credible estimate (at 
least in theory) of the catchment-mean rates. 

The catchment LSM uses three non-traditional bulk moisture variables to 
describe the catchment moisture state, one representing equilibrium condi- 
tions associated with the water table distribution, and the others represent- 
ing non-equilibrium conditions near the surface. The LSM currently uses 
the TOPMODEL formulation of Beven and Kirkby [ 1979 ] to establish the 
water table distribution. The limitations of TOPMODEL in low topography 
regions are diminished through the use of the non-equilibrium moisture vari- 
able in low topography regions, the model reduces to a more traditional 
layer-type LSM. 

The catchment strategy has several advantages. The two most important 
are: (i) the explicit treatment of subgrid soil moisture variability allows for 
more realistic formulations of evaporation and (especially) runoff; and (n) the 
use of the catchment as the land surface element allows a strong connection 
between LSM structure and established hydrological models, which will be 
helpful in the future evolution of our LSM. The overall strategy, which builds 
on the work of others who have recognized the importance of the runoff 
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formulation in LSMs [Famiglietti and Wood , 1991; Liang et al. , 1994; Stieglitz 
et al., 1997], avoids many of the problems limiting standard representations 
of land surface processes. 

Although the catchment strategy’s conceptual advantages are important 
in themselves, the strength of the approach is best demonstrated through 
an application of the fully coded model in a real setting, with observed me- 
teorological variables for forcing the model and observed water fluxes for 
evaluating the model’s performance. Such an analysis is presented in Part 2 
[Ducharne et al., this issue]. 
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Appendix A: Equations Used to Partition the 
Catchment 

In the discussion below, the average equilibrium root zone moisture at states 
(A), (B), and (C) in Figure 4 will be referred to as 0^, 0^, and 9^, the 
catchment deficits at states (A), (B), (C) and (D) will be referred to as 
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M d , M (B) , M (C) , and M (D) , and the equilibrium pdfs at states (A), (B) and 
(C) will be referred to as /($), /(b)(0) and /(c)(0)- The average root zone 
moisture at state (D) is, by definition, 0 wi i t . 

Also, a dimensionless root zone excess, 0 exc , is used below. It is computed 
by dividing the root zone excess prognostic variable, M TZ , by the water hold- 
ing capacity of the root zone. The mean root zone wetness at a given time, 
corrected for the root zone excess, is thus 0(A) + 0 exc . 

The equations for A sat , A tr , and ,4 wilt vary according to where M D and 0 exc 
position the non-equilibrium pdf of soil moisture in Figure 4. If M D < 
so that a water table is active, and if the root zone excess does not bring any 
part of the pdf below the wilting point, then we compute 

/ OO 

Z(0 - 0 exc )d0 

Atr = 1 • A sat 


A w iit = 0. 


If Md < M but the root zone excess does bring part of the pdf below the 
wilting point, we assume the wilting fraction to vary linearly between 0 (at 
0(a) + 0 exc = 0(c)) and 1 (at 0(a) + 0 exc = 0 w n t ): 


a _ 0(C) ~ 0 ( A ) - 0 exc 

^ wilt — r — a 

tqc) — wilt 

/ OO 

f(0 - 0 exc )d9 


Atr — I- A sat A w itt- 
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If M D > then no water table is active, and the “equilibrium” 

average root zone moisture associated with Mp is computed assuming a 
linear variation of 0 ^ with M D : 

-5— _ Q , (s— Q \ - M d 

0(A) - 0 wilt + (0(B) - 0wilt) _ m(b) ■ (15) 

Adding 9 exc to 9 (a) in this “no water table” case will bring the average root 
zone moisture into one of three regimes. First, if 9^ + 9 exc > 9^ } (i. e ., the 
average root zone moisture is above 0 (b)), the areas are computed with 

^ Sat = Jl + ^( B ) ~ ^(A) — 0exc)d9 

-^tr = 1 • -^-sat 

-<4 wilt = 0. ( 16 ) 

Second, if 0 (c) < 0(A) +0 e *c < 0 (b), then the final pdf lies between states (B) 
and (C), and the areas are calculated with 

A - a(Q , ^(A) + 0exc - 0(C) ( ,,(B) ,fC)\ 

^sat - A sat + 1 ' Msat ~ A ( J) 

y (B) ~ 0(C) V ' 

Atr = 1. -4sat 

•^wilt = 0, (17) 

where 

AS? = /(B)(0)d0 (18) 

and 

AS = f™ f( C )(9)d9. (19) 
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Finally, if 0<a) 4- 0 ex c < 0(C), wilting occurs over a portion of the catchment. 
We compute: 

a _ 0(C) “ 0(A) - 0exc 
^wilt — q a 

0(C) ~ tf w ;it 

A S at = (1 — -^wilt) 

Atr = 1 A sa t (20) 
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Figure Captions 


Fig. 1 Illustration of catchment delineation in southwest North America. 
Overlain on the plot is a l°xl° grid and a single 4°x5° grid cell. 

Fig. 2 Illustration of the calculation of “catchment deficit”, a bulk water 
variable used by the new LSM. 

Fig. 3 Illustration of “root zone excess” at a point. 

Fig. 4 Probability density functions (pdfs) of equilibrium root zone soil mois- 
ture at various catchment deficits in the catchment; see text for details. 

Fig. 5 Thermal heat flux calculations in the catchment LSM. The same areal 
fraction of snowcover is assumed to lie over the saturated fraction (S), 
the transpiration fraction (T), and the wilting fraction (W). In each 
fraction, the ground surface exchanges heat with the snowcover and 
with the overlying atmosphere separately, and the snow surface itself 
exchanges heat with the atmophere. The three surface fractions ex- 
change heat separately with an underlying ground layer, which in turn 
exchanges heat with lower ground layers. 
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Figure 1: Illustration of catchment delineation in southwest North America. 

Overlain on the plot is a l°xl° grid and a single 4°x5° grid cell. 
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Moisture profile at an arbitrary point in the catchment: 
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Moisture profiles across catchment: 
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Figure 2: Illustration of the calculation of “catchment deficit”, a bulk water 
variable used by the new LSM. 
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2ND SOIL MOISTURE PROGNOSTIC VARIABLE: 
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Figure 3: Illustration of “root zone excess” at a point. 
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(D): complete wilting, 
(delta function) 


(B): mean water table 



Figure 4: Probability density functions (pdfs) of equilibrium root zone soil 
moisture at various catchment deficits in the catchment; see text for details. 
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Figure 5: Thermal heat flux calculations in the catchment LSM. The same 
areal fraction of snowcover is assumed to lie over the saturated fraction (S), 
the transpiration fraction (T), and the wilting fraction (W). In each fraction, 
the ground surface exchanges heat with the snowcover and with the overlying 
atmosphere separately, and the snow surface itself exchanges heat with the 
atmophere. The three surface fractions exchange heat separately with an 
underlying ground layer, which in turn exchanges heat with lower ground 
layers. 
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Popular Summary: 


“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. Part 
1: Model Structure”, by R. Koster, M. Suarez, A. Ducharne, M. Stieglitz, and P. 
Kumar 

“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. Part 
2: Parameter Estimation and Model Demonstration”, by A. Ducharne, R. Koster, M. 
Suarez, M. Stieglitz, and P. Kumar 

In the standard approach to modeling land surface processes in general circulations 
models (GCMs), unresolved heterogeneity in soil moisture is ignored. This severely 
compromises the accuracy of simulated large-scale surface fluxes, since many critical 
hydrological processes respond nonlinearly to spatial variations in soil moisture. 

In our two-part paper, we present a new land surface model designed to overcome 
this deficiency. Part 1 provides the overall modeling strategy and details in the model 
formulation. The strategy first involves using highly resolved topographical data 
to separate the continental surface into a mosaic of hydrological catchments. The 
irregular catchment grid is quite distinct from the quasi-rectangular grid employed 
by most land surface models (LSMs). Simple transformations convert the forcing 
from the overlying regular atmospheric grid to the irregular catchment grid. 

The strategy then involves the diagnosis, for each catchment, of a soil moisture 
distribution function from bulk measures of catchment moisture content combined 
with the statistics of the catchment’s topography. The distribution function is used to 
separate the catchment into three hydrologically-distinct areal fractions: the fraction 
that is saturated at the surface, the fraction that is subsaturated but still transpiring 
freely, and the fraction this is undergoing wilting. Evaporation and runoff processes 
are calculated independently over each areal fraction, using formulations relevant to 
each fraction’s specific hydrological regime. 

Part 2 addresses the appropriateness of the simple functional forms used to ap- 
proximate the relevant complex theoretical equations, and it demonstrates the fea- 
sibility of establishing model parameter values on a global scale. More importantly, 
it demonstrates that the model described in Part 1 does indeed work well ^ with a 
minimal amount of calibration, it can reproduce observed large-scale evaporation and 
runoff rates in the Red-Arkansas Basin. Additional applications of the model show 
that topographical variations, as modeled, do have a significant impact on simulated 
evaporation and runoff rates. 
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The most important advantage of the overall approach is the potential for im- 
proved surface fluxes through the explicit consideration of soil moisture heterogene- 
ity. The strategy boasts additional advantages, though, including: (a) the potential 
for the straightforward application of river routing schemes; (b) the diagnosis of soil 
moisture contents that are more amenable to validation with observations; and (c) 
the potential for topography-based disaggregation of GCM forcing. In addition, the 
specification of the catchment as the fundamental land surface element provides a 
more direct link to basic hydrological science, for which the catchment is a natural 
unit of study. 
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“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. 
Part 1: Model Structure”, by R. Koster, M. Suarez, A. Ducharne, M. Stieglitz, 
and P. Kumar 

“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. 
Part 2: Parameter Estimation and Model Demonstration”, by A. Ducharne, R. 
Koster, M. Suarez, M. Stieglitz, and P. Kumar 

\ 

\ 

Statement of Significance 

QUESTION: Most land surface models (LSMs) used with general circulation 
models (GCMs) today do not include adequate treatments of soil moisture’s 
spatial variability, and this severely limits their production of realistic surface 
fluxes. How can this fundamental deficiency be addressed? 

APPROACH: Our two papers describe the framework and application of a new 
kind of land surface parameterization. The hydrological catchment is 'defined 
as the basic land surface element, and within each catchment, the spatial het- 
erogeneity of soil moisture is diagnosed from bulk moisture variables and the 
catchment’s topography. This spatial heterogeneity is then employed in a more 
sensible treatment of subsurface moisture transfer, runoff generation and evap- 
oration production. 

SIGNIFICANCE AND IMPLICATIONS OF FINDINGS: The explicit consid- 
eration of subgrid-scale soil moisture heterogeneity and its effects on the surface 
energy and water budgets, as described in Part 1, should in principle lead to an 
improved simulation of land surface processes and thus to a more reliable simu- 
lation of climate and its sensitivity. The applications of the model described in 
Part 2 demonstrate that the model can successfully reproduce observed fluxes 
and is computationally viable for global application. 


RELATIONSHIP TO MTPE SCIENCE PLAN 

The improved simulation of land surface processes should benefit numerical 
climate studies. The strategy described in the paper thus has direct relevance to 
MTPE emphases on (1) seasonal-to-interannual climate prediction, (2) changes 
in long term climate, and (3) landcover and land use change. 
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“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. 
Part 1: Model Structure”, by R. Koster, M. Suarez, A. Duchame, M. Stieglitz, 
and P. Kumar 

“A Catchment-Based Approach to Modeling Land Surface Processes in a GCM. 
Part 2: Parameter Estimation and Model Demonstration”, by A. Ducharne, R. 
Koster, M. Suarez, M. Stieglitz, and P. Kumar 


Statement of Significance 
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fluxes. How can this fundamental deficiency be addressed? 

APPROACH: Our two papers describe the framework and application of a new 
kind of land surface parameterization. The hydrological catchment is defined 
as the basic land surface element, and within each catchment, the spatial het- 
erogeneity of soil moisture is diagnosed from bulk moisture variables and the 
catchment’s topography. This spatial heterogeneity is then employed in a more 
sensible treatment of subsurface moisture transfer, runoff generation and evap- 
oration production. 

SIGNIFICANCE AND IMPLICATIONS OF FINDINGS: The explicit consid- 
eration of subgrid-scale soil moisture heterogeneity and its effects on the surface 
energy and water budgets, as described in Part 1, should in principle lead to an 
improved simulation of land surface processes and thus to a more reliable simu- 
lation of climate and its sensitivity. The applications of the model described in 
Part 2 demonstrate that the model can successfully reproduce observed fluxes 
and is computationally viable for global application. 


RELATIONSHIP TO MTPE SCIENCE PLAN 

The improved simulation of land surface processes should benefit numerical 
climate studies. The strategy described in the paper thus has direct relevance to 
MTPE emphases on (1) seasonal-to-interannual climate prediction, (2) changes 
in long term climate, and (3) landcover and land use change. 
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